Accurately model the Kuramoto--Sivashinsky dynamics with holistic discretisation

We analyse the nonlinear Kuramoto--Sivashinsky equation to develop accurate discretisations modeling its dynamics on coarse grids. The analysis is based upon centre manifold theory so we are assured that the discretisation accurately models the dynamics and may be constructed systematically. The theory is applied after dividing the physical domain into small elements by introducing isolating internal boundaries which are later removed. Comprehensive numerical solutions and simulations show that the holistic discretisations excellently reproduce the steady states and the dynamics of the Kuramoto--Sivashinsky equation. The Kuramoto--Sivashinsky equation is used as an example to show how holistic discretisation may be successfully applied to fourth order, nonlinear, spatio-temporal dynamical systems. This novel centre manifold approach is holistic in the sense that it treats the dynamical equations as a whole, not just as the sum of separate terms.Comment: Without figures. See http://www.sci.usq.edu.au/staff/aroberts/ksdoc.pdf to download a version with the figure

Use the information dimension, not the Hausdorff

Multi-fractal patterns occur widely in nature. In developing new algorithms to determine multi-fractal spectra of experimental data I am lead to the conclusion that generalised dimensions $D_q$ of order $q\leq0$, including the Hausdorff dimension, are effectively \emph{irrelevant}. The reason is that these dimensions are extraordinarily sensitive to regions of low density in the multi-fractal data. Instead, one should concentrate attention on generalised dimensions $D_q$ for $q\geq 1$, and of these the information dimension $D_1$ seems the most robustly estimated from a finite amount of data.Comment: 11 page

Computer algebra models the inertial dynamics of a thin film flow of power law fluids and other non-Newtonian fluids

Consider the evolution of a thin layer of non-Newtonian fluid. I model the case of a nonlinear viscosity that depends only upon the shear-rate; power law fluids are an important example, but the analysis is for general nonlinear dependence upon the shear-rate. The modelling allows for large changes in film thickness provided the changes occur over a large enough lateral length scale. The modelling is based on two macroscopic modes by fudging the spectrum: here fiddle the surface boundary condition for tangential stress so that, as well as a mode representing conservation of fluid, the lateral shear flow u ∝ y is a neutral critical mode. Thus the resultant model describes the dynamics of gravity currents of non-Newtonian fluids when their flow is not very slow. For an introduction I first report on an analogous case of nonlinear diffusive dissipation